A measurement of the broad-band spectrum of XTE Jll 18+480 
with BeppoSAX and its astrophysical implications 
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ABSTRACT 

We report on results of a target of opportunity observation of the X-ray transient XTE 
J1118+480 performed on 2000 April 14-15 with the Narrow Field Instruments (0.1-200 keV) of 
the BeppoSAX satellite. The measured spectrum is a power law with a photon index of ~ 1.7 
modified by an ultrasoft X-ray excess and a high-energy cutoff above ~ 100 keV. The soft excess 
is consistent with a blackbody with temperature of ~ 40 eV and a low flux, while the cut-off 
power law is well fitted by thermal Comptonization in a plasma with an electron temperature of 
~ 10 2 keV and an optical depth of order of unity. Consistent with the weakness of the blackbody, 
Compton reflection is weak. Though the data are consistent with various geometries of the hot 
and cold phases of the accreting gas, we conclude that a hot accretion disk is the most plausible 
model. The Eddington ratio implied by recent estimates of the mass and the distance is ~ 10~ 3 , 
which may indicate that advection is probably not the dominant cooling mechanism. We finally 
suggest that the reflecting medium has a low metallicity, consistent with location of the system 
in the halo. 
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1. Introduction 

The high Galactic latitude source (I = 157.7°, 
b = 62.3°) XTE Jll 18+480 was discovered in 
the 2-10 keV energy band with the All-Sky Mon- 
itor (ASM) aboard the Rossi X-Ray Time Ex- 
plorer (RXTE) satellite on 2000 March 29 as a 
weak (~ 40 mCrab) X-ray source. Retrospec- 
tive analysis showed that its flux was slowly ris- 
ing since 2000 March 5 (Remillard et al. 2000) and 
that another outburst had occurred in 2000 Jan- 
uary 2-29, with a similar peak flux on January 6 
(Remillard et al. 2000), see Figure 1. The source 
was also detected in the 20-100 keV band with 
the CGR O/BATSE experiment (Wilson & McCol- 
lough 2000): the flux was ~ 60 mCrab on March 
26 and ~ 110 mCrab on January 11. The energy 
spectra derived from the ASM and BATSE data 
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were described by a power law with photon index 
of <~ 1.8 and ~ 2.1, respectively. These power law 
slopes are typical of black hole binaries (e.g., Cyg 
X-l) in the hard state. 

A radio counterpart of XTE J1118+480 was 
discovered by Pooley & Waldram (2000) on March 
31 with a flux density of 6.2 mJy at a frequency of 
15 GHz. Also, the source was detected in the op- 
tical and near-infrared bands (Garcia et al. 2000; 
Chaty et al. 2000; Wren & McKay 2000) with 
magnitudes of 12.12, 11.75 and 11.06 in the J, 
H, K bands, respectively, on 2000 April 4.2 (UT), 
and of 13.11 in the V band on March 28.206. 
Optical spectroscopy performed on March 31.0 
shows Ha emission with an FWHM (2200 km 
s _1 ) that is approximately equal to the largest 
found in black hole X-ray novae during qui- 
escence (Garcia et al. 2000). Observations per- 
formed with the Hubble Space Telescope (Haswell, 
Hynes, & King 2000a) on 2000 April 8 show an 
Ha absorption line with full width corresponding 
to radial velocities likely > 10 4 km s _1 , which are 
suggestive of a massive accretor. 

The optical counterpart of XTE Jl 118+480 is 
found at RA. (2000) = ll h 18 m 10 s .85 and Decl. 
(2000) = +48°02'12".9 (Uemura et al. 2000). 
A star with a 18.8 R magnitude in the USNO 
A1.0 and A2.0 catalogues is within 2" of that di- 
rection, and it has been proposed as the quies- 
cent optical counterpart of the source (Uemura 
et al. 2000). Its low magnitude is consistent 
with the LMXB nature of the transient source 
(Tanaka & Shibazaki 1996). An interesting prop- 
erty of the system in outburst is the very low 
X-ray-to-optical flux ratio (~ 5) with respect 
to the typical value of ~ 500 found in LMXBs 
(van Paradijs & McClintock 1995). 

Cook et al. (2000) reported the discovery of a 
sinusoidal periodic variation of the V-band light 
curve with a period of 0.1706 days, while a value 
of 0.17078 ± 0.00004 days has been obtained by 
Uemura et al. (2000) . No evidence of periodicity 
was found in the soft X-rays (Hynes et al. 2000). 
The above optical modulations have been inter- 
preted as a superhump (Kuulkers 2001), while the 
true binary period (P = 0.1703 days) has been 
determined by McClintock et al. (2000, 2001) 
and Wagner et al. (2000) through radial velocity 
studies of the optical counterpart in quiescence. 
Rcvnivtsev, Sunyaev, & Borozdin (2000) report 



the discovery of a QPO feature at 0.085 Hz in the 
X-ray light curves from the Proportional Counter 
Array (PCA) aboard RXTE, while the evolution of 
this QPO for a long time has been investigated by 
Wood et al. (2000). Besides the QPO, the Power 
Spectral Density (PSD) shows a strong band- 
limited noise component (40% fractional rms), 
also typical of BHCs in their low state. No signifi- 
cant variability was detected at higher frequencies. 
Similar variations and QPOs were observed in the 
optical/UV band (Haswell et al. 2000b). 

The definite confirmation of the black hole na- 
ture of the source has been provided by mea- 
surements of the mass function being as large as 
~ 6M Q (McClintock et al. 2000, 2001a, hereafter 
MOla; Wagner et al. 2000), rivaled only by that of 
V404 Cyg. In §2.1 below, we discuss the black hole 
mass and inclination implied by those and other 
measurements. 

Given the broad energy band of the instru- 
mentation on board BeppoSAX (0.1-200 keV) and 
the high Galactic latitude of the source, XTE 
J1118+480 could be searched for an ultrasoft com- 
ponent during a hard spectral state with unprece- 
dented sensitivity. As soon as the source be- 
came observable with BeppoSAX , we started a 
Target Of Opportunity (TOO) observational pro- 
gram devoted to study the spectral properties of 
the source. Here we present results of the first ob- 
servation performed on 2000 April 14-15. In figure 
1 we show the epochs of the source observations 
with BeppoSAX. 

2. Parameters of the System 

2.1. Masses, inclination and distance 

McClintock et al. (2000), MOla and Wagner et 
al. (2000) have recently detected absorption lines 
of an M0-M1 V star in their red optical spectra 
obtained when the object was in a quiescent state. 
The weighted average of their measurements of the 
velocity amplitude for the secondary is Ki — 697± 
9 km s _1 , which yields the mass function of 



f(M 2 ) = 



(M x sini) 



6.0 ± 0.2 M 



© 



(1) 



for P — 0.17 days (where % is the inclination). 

At the same time, Garnavich (2000; see also 
MOla) found a double sinusoidal variation with a 
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Fig. 1 . — One-day-average count rates from XTE 
Jl 118+480 detected with the RXTE/ ASM during 
the outburst. (From the Internet public archive at 
xtc.mit.cdu/XTE/asmlc/ASM.html.) The Beppo- 
SAX observations are marked with the arrows. 
The observation analyzed here is the earliest one. 

peak-to-peak amplitude of 0.15 mag in the /-band 
light curve phased to the 0.17-day period. This 
suggests that the secondary is tidally distorted. 

In systems with mass ratios of q = M-^fM-i > 
1.25, the average density of a Roche lobe-filling 
secondary is determined solely by the orbital pe- 
riod, which implies simple period-mass-radius- 
spectral type relations for lower main sequence 
(MS) secondaries (Frank, King, & Raine 1992). 
In the case of XTE 1118+480, the orbital pe- 
riod implies a spectral type of M2 and a mass 
of ~ 0.45 M© for a MS secondary with near-solar 
metallicity, whereas a star with lower metallicity 
would have somewhat earlier spectral type (e.g. 
Beuermann ct al. 1998). An evolved secondary 
would have a spectral type significantly later than 
that of a zero-age MS star with the same density 
(i.e., in a binary with the same orbital period) and 
a lower mass (e.g. Baraffe & Kolb 2000). If the 
secondary is a MS star with M 2 ~ 0.45 M Q , equa- 
tion (1) yields q > 15, and M x sin 3 i > 6.8 M Q . 
The lack of dips or eclipses for a Roche-lobe filling 
secondary, in particular in the X-ray light curve 
(Revnivtsev et al. 2000), yields an upper limit for 
the inclination of i < 80°, and M x > 7.1 M Q . 

The light curve presented in fig. 1 of MOla 
shows some deviations from a purely ellipsoidal 
light curve. First, it shows some extra light at 



~ 0.05P after the inferior conjunction of the sec- 
ondary, which, according to MOla can be ex- 
plained by a bright spot at the edge of accretion 
disk. There is also a phase lag of 0.05 relative to 
the spectroscopic ephemeris, which makes the el- 
lipsoidal nature of the light curve somewhat prob- 
lematic. On the other hand, MOla found a sig- 
nificant contamination of the visual light by an 
accretion disk, which likely contaminates also the 
infrared light. In this case, eclipse effects in an 
extended disk can possibly shift the phase of the 
light curve. 

Thus, we follow MOla and assume that the light 
curve is dominated by ellipsoidal effects. The am- 
plitude of the double sinusoid is then positively 
correlated with q, i, and the Roche-lobe filling fac- 
tor. Moreover, the two minima of the cllipticity 
effect arc in general unequal, the one correspond- 
ing to the superior conjuntion of the distorted star 
being deeper than the other. The magnitude dif- 
ference between the two minima is a function of i 
and the filling factor. In particular, for a Roche- 
lobe filling star, this difference should be about 0.1 
mag for i > 80° decreasing almost linearly to at 
i~40°. 

The light curve presented by MOla shows the 
shallower minimum at the inferior conjunction (as 
expected) of Al ~ 0.11, and the magnitude differ- 
ence between the two minima of ~ 0.07, indicat- 
ing i > 60°. On the other hand, for g > 15 and 
a Roche-lobe filling star with convective envelope, 
we estimate the amplitude of the shallower mini- 
mum to be Al > 0.35 sin i mag, which suggests 
much lower inclination, i ~ 34°. This shows that 
the amplitude of the ellipsoidal variability is re- 
duced due to contamination of the 7-band flux by 
some additional contribution from other source (s), 
e.g. an accretion disk. 

MOla found that the secondary contributes 42% 
(K5V) to 26% (M1V) of the total flux at 5900 A. 
To evaluate a contribution of the accretion disk at 
the /-band, we assume that the local disk spec- 
trum corresponds to either F v or F\ being con- 
stant. Assuming an M1-M2 secondary, we esti- 
mate the contribution of the disk as ~ 30% and 
~ 60%, in the first and second case, respectively. 

Then, the following equations can be used to 
determine q, i, and the binary separation, a. Ke- 



pier's third law reads, 

M x + M 2 = 0.01343 . |^ M . (2) 
(P/l day) 2 

The secondary's radial velocity amplitude, K 2 = 
697 km/s, corresponds to 

asini = 2.34^^R . (3) 

q 

The full (peak-to-peak) amplitude produced by 
the tidal distortion is 

A/ = 4.3<7(— J sin 2 i mag, (4) 

where i? t corresponds to the "transverse radius" 
of the secondary, which is about 96% of the 
mean radius, R 2l of its Roche lobe. The nu- 
merical coefficient in equation (4) was evaluated 
from synthetic /-band light curves generated with 
the Wilson-Devinney code (Wilson 1990) using 
a gravity-darkening exponent of 0.32, a linear 
limb-darkening coefficient of 0.6, and using the 
shallower minimum in the synthetic light curve. 
Assuming a MS secondary with M 2 ~ 0.45M Q 
and R 2 — O.45R0, we obtain a — 3.5R0, q = 44, 
i = 43°, M x = 20M Q , and a = 2.8R Q , q = 22, 
i = 61°, Mx = 10M© in the cases with con- 
stant F v and F\, respectively. The first solution 
produces equal minima, which is not observed, 
whereas the second one results in a too shallow 
second minimum. 

On the other hand, the presence of extra light 
near the first minimum will reduce its depth, 
whereas the second minimum may be not affected. 
Using only the depth of the second minimum, we 
get a = 3.13R , q = 31, i = 50°, M x = 14M Q 
in the first case, and a solution similar to that of 
MOla, a = 2.53R Q , q = 16, i = 80°, M x = 7.1M , 
in the second case. Note that the first solution still 
produces equal minima in the light curve. 

We note that the determinations of the binary 
parameters by MOla and above supersede analo- 
gous constraints of Dubus et al. (2001), obtained 
before the recent measurements of K 2 . Those 
authors used instead the peak-to-peak separation 
of very broad double-peaked Hi Balmer and Hen 
4686 A emission lines, measured by them as At; w 
1270, 1470 km/s, respectively. They then assumed 
that Aw represents the projected Keplerian veloc- 
ity of the outer disk, and used a theoretical value 



of its radius, i?d- For a Keplerian disk with a 
power-law emissivity, the relation, 



Q~' = 



Av 
2K~ 2 



Rd 
a 



(5) 



is satisfied (e.g. Smak 1981). However, with the 
measured K 2l 1 + q^ 1 w 0.8i?d/ a , 1.1-ffd/a for the 
Balmer and Hen lines, respectively. This implies 
Rd ~ a for any value of q, which then appears to 
invalidate their adopted assumption of a Keplerian 
outer disk. 

Summarizing, the present most likely range of 
the mass and inclination is 

60°<i<80°, 1O<M X /M <7. (6) 

The lower and upper limit on i follow from the un- 
equal light-curve minima and the lack of eclipses, 
respectively. 

The presence of the Ml-2 V star with V ~ 20.4 
(MOla) implies the distance of 



d~ (1.3-1.9) kpc, 



(7) 



where the error also takes into account the un- 
certainty in the interstellar absorption (see be- 
low). We also assumed My ~ 8.96-9.9 for such 
a MS star with near-solar metallicity (Baraffe & 
Chabrier 1996). Note that d will become lower for 
lower metallicity (Beuermann et al. 1998). 

2.2. The hydrogen column density 

An accurate determination of the hydrogen col- 
umn density, Ah, is of key importance for deter- 
mining the correct spectral model for the source 
emission. Given the high Galactic latitude of the 
source, the total Galactic column density along 
the source direction is likely to be close to Ah, 
as well as it represents an upper limit to it. The 
Galactic column density of HI has been measured 
in radio in 5 directions within 1° of the direction 
of the optical counterpart, yielding the range of 
A H , G = (1.28-1.44) x 10 20 cm" 2 (Dickey & Lock- 
man 1990, as implemented in the ftools soft- 
ware package). The measurement in the direc- 
tion closest (0.29°) to that of XTE J1118+480 
gives 1.28 x 10 20 cm -2 , and the average weighted 
by the angular distance from the source direc- 
tion is 1.34 x 10 20 cm -2 . We note that, at these 
low values of Ah, absorption is effective only at 



soft X-rays, at which Hi is indeed the main ab- 
sorbing species. A lower Galactic column den- 
sity (0.75 x 10 20 cm~ 2 ) along the line of sight 
is derived following Schlegel, Finkbeiner & Davis 
(1998), who obtained full-sky dust maps on the 
basis of the COBE and IRAS satellite data. The 
uncertainty in their reddening maps is estimated 
to be about 16%. Taking into account that the Ah 
value along the line of sight, obtained by interpo- 
lation from these and Hi maps, can differ from the 
real value by a factor ~ 2 (Faison et al. 1998), we 
can state that the above Ah estimates are consis- 
tent each other. 

Then, Dubus et al. (2001) have found 3 weak 
absorption features (Can 3933 A doublet and 
3968 A), interpreted as coming from 3 interstel- 
lar clouds in the line of sight moving at different 
velocities. From best-fitting, they obtain the to- 
tal Acaii — 1.4 x 10 12 cm~ 2 (without giving its 
uncertainty). They then convert this value to Ah 
using an average range of (log(AH/Ac a ii)) ~ 8.1- 
8.5 obtained for high-altitude lines of sight (Sem- 
bach & Danks 1994). This yields Ah in the range 
of (1.76-4.4) x 10 20 cm -2 . We note, however, that 
it is sufficient that the actual Nn/Ncaii ratio in 
the direction of XTE Jl 118+480 is just - 25% 
less than the lower limit of the above average ra- 
tio range to achieve the consistency of Ah from 
radio mapping of Hi with that from Can. In fact, 
the H and K lines of Can correlate poorly with 
-/Vh in general and are thus considered not suitable 
for the purpose of accurate determination of .Ah 
(Dickey & Lockman 1990 and references therein). 

On the other hand, Hynes et al. (2000) obtained 
A H in the range of (0.35-1.15) x 10 20 cuT 2 by fit- 
ting their 0.1-0.2 keV EUVE spectrum by a power 
law with a range of indices. We note that such 
narrow-band fits appear to be easily subject to 
systematic errors due to the possibility that the 
actual spectrum is not a power law. Our Beppo- 
SAX data include this range, and we have found 
Ah ~ (0.7-1.5) x 10 20 cm~ 2 by fitting a range of 
models. Power-law models gave lower values of 
Ah, consistent with Hynes et al. (2000), but mod- 
els with blackbody and Comptonized blackbody 
spectra gave higher Ah, consistent with the mea- 
surements of Dickey & Lockman (1990). Given the 
results above, we consider the most likely range of 
A H to be ~ 1-1.5 x 10 20 cm~ 2 . 



3. The Observation and Spectral Data 

The observation was performed with the Beppo- 
SAX Narrow Field Instruments (NFIs). The 
RXTE/ ASM light curve of XTE Jl 118+480 in 
the 2-10 keV energy band with the mark of the 
BeppoSAX observations is shown in Figure 1. The 
NFIs embody a Low Energy Concentrator Spec- 
trometer (LECS, 0.1-10 keV, Parmar et al. 1997), 
2 Medium Energy Concentrators Spectrometers 
(MECS, 1.3-10 keV, Boella et al. 1997), a High 
Pressure Gas Scintillator Proportional Counter 
(HPGSPC, 3-100 keV, Manzo et al. 1997), and a 
Phoswich Detection System (PDS, 15-300 keV, 
Frontera et al. 1997). The LECS and MECS 
have imaging capabilities, while the HPGSPC and 
PDS are collimated detectors with a Field of View 
(FOV) of 1° and 1.3°, respectively. The PDS 
instrument makes use of rocking collimators for 
background monitoring. 

We observed the source from 2000 April 
14.489 UT to 15.528 UT. The source exposure 
times were 21197 s, 29758 s, 42138 s, and 20146 s 
for the LECS, MECS, HPGSPC and PDS, respec- 
tively. Useful data were selected from time inter- 
vals that met the following criteria: the satellite 
outside the South Atlantic Geomagnetic Anomaly, 
the elevation angle above the Earth limb by > 5°, 
dark Earth (for LECS), stabilized high voltage 
supplies. 

The source was detected with all the NFIs, from 
-0.1 to 200 keV. The LECS and MECS source 
spectra were extracted from a region with a ra- 
dius of 8' around the centroid of the source image. 
As background spectra we used standard files ob- 
tained from the observation of blank fields. The 2 
MECS spectra were equalized and co-added. We 
used background spectra accumulated from dark 
Earth data for the HPGSPC (note that the instru- 
ment's collimator was kept on the source for the 
entire observation). The background level of the 
PDS was estimated by swapping its collimators 
off source every 96 s. The energy bands used for 
spectral fitting were limited to those where the re- 
sponse functions were best known, i.e, 0.12-4 keV, 
1.7-10 keV, 7-29 keV, and 15-200 keV, for the 
LECS, MECS, HPGSPC and PDS, respectively. 
The count rate spectra were analyzed using the 
xspec software package (Arnaud 1996). A sys- 
tematic error of 1% was added in quadrature to 
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Fig. 2.— Top panel: The 0.1-200 kcV count 
rate spectrum of XTE J1118+480. Bottom panel: 
Ratio between the count rate spectrum and the 
power-law model absorbed by the Galactic TVh- 
A soft X-ray excess and a high-energy cutoff are 
clearly highly significant. 



the statistical uncertainties of the spectral data, 
on the basis of the calibration results obtained 
with the Crab Nebula, that was observed 4 days 
before our observation (10 April 2000). We al- 
lowed for free normalization of the instruments in 
multi-instrument fits with respect to MECS. For 
clarity of display, the unfolded spectra from multi- 
instrument fits were renormalized to the MECS 
level. 

The results below have been obtained with the 
spectra averaged over the entire observation. We 
have also performed a spectral analysis on differ- 
ent time intervals, but have not found any statis- 
tically significant variation of the best-fit parame- 
ters with time. 

The quoted errors for the spectral parameters 
correspond to 90% confidence for one parameter 
(A% 2 = 2.71). The elemental abundances are with 
respect to those of Anders & Ebihara (1982), and 
the opacities are from Morrison & McCammon 
(1983). In calculations of Compton reflection and 
in Comptonization models in planar geometry, we 
assume an inclination of i — 70° (see §2.1). 

4. Results 

The 0.1-200 keV count rate spectrum obtained 
with the NFIs (see Figure 2, top panel) was fit with 
different models. A simple power-law model (pl) 



unacceptable fit with the photon index, T = 1.81, 
and a very large x 2 l v — 1143/231, caused mostly 
by strong systematic residuals at the lowest and 
highest energies, corresponding to a soft X-ray 
excess at < 0.5 keV and a high-energy cutoff at 
> 100 keV, respectively, as shown in the bottom 
panel of Figure 2. The high energy cutoff is also 
apparent from figure 3, where we show the ratio 
of the XTE J1118+480 spectrum measured by the 
PDS with the corresponding one of the Crab Neb- 
ula, which was observed for cross-calibration pur- 
poses 4 days before our observation. Replacing 
the power law with an e-folded power law (cut- 
OFFPl) and maintaining the iVn value fixed to 
1.28 x 10 20 cm~ 2 only marginally improve the fit 
(x 2 /^ = 1136/230): the systematic residuals at 
low and high energies are still there (see Figure 2), 
and the best fit cut-off energy is very high and un- 
constrained. Only allowing a free JVh the fit does 
improve significantly, to x 2 l v = 319/229, but only 



for an TVjj = 0.66 



+0.03 
-0.03 



10 20 cm 2 , about a half 



with JVi 



H 



1.28 x 10 20 cm" 2 (see §2.2) gives an 



of Nn,G, even if consistent with the TVh based on 
the dust IR maps (§2.2). However still there re- 
main strong systematic residuals at the lowest and 
highest energies, as confirmed by the still unac- 
ceptable x 2 ■ Thus we added to the previous model 
a blackbody spectrum (bb). With this additional 
component the fit quality significantly increased 
{x 2 /v = 272/227), with a chance probability that 
the lower x 2 is due to chance of 1.4 x 10~ 8 . Even 
if the CUTOFFPL model component is not the best 
description of the high energy cut-off (see below), 
it represents a great improvement with respect to a 
PL component, that gives a much worse and unac- 
ceptable quality of the fit { X 2 /v = 388/228). This 
result confirms the statistical significance of the 
high energy cut-off in our data. The best fit model 
with BB plus CUTOFFPL provides an iVn close to 
./Vh.g with a rather low temperature {kTy,^ ~ 50 
cV) of the BB spectrum, see Table 1. Figure 4 
shows the best fit curve along with the residuals 
of the data to the model. 

Still, even the last fit shows strong systematic 
residuals at the highest energies, as an e- folded 
power law does not reproduce the actual shape 
of the high-energy cutoff in the data. Therefore, 
we have tested models with the main continuum 
component given by Comptonization, a common 
physical process in compact sources (see Zdziarski 
2000 for a review). We first consider the thermal- 



Table 1 
Fit results of the 0.1-200 keV spectrum of XTE J1118+480 
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Fig. 3. — Ratio between the count spectrum of 
XTE Jl 118+480 and that of Crab observed with 
BeppoSAX 4 days before our TOO measurement. 
The break in the ratio at ~ 100 keV is apparent. 
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Fig. 4.— Top panel: The 0.1-200 keV count rate 
spectrum of XTE J1118+480, superposed with the 
best fit model consisting of an absorbed BB plus 
CUTOFFPL (see text). Bottom panel: Ratio be- 
tween the count spectrum and the model. 



Compton model (comptt) of Titarchuk (1994), 
in which a Wien spectrum undergoes scattering 
in either spherical or slab geometry, and rela- 
tivistic effects are taken into account. We have 
found that fits with very similar electron temper- 
ature, kT, and \ 2 are obtained for either geom- 
etry. Similarly to the case of an e-folded power 
law, the absorbed comptt model yields an im- 
plausibly low N H (0.51±g;?| x 10 20 cm" 2 ) with a 
X /v = 283/228, similar to that found with the 
last model (bb plus CUTOFFPL). Adding an ab- 
sorbed bb component, with Tbb equal to the tem- 
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Fig. 5.— Top panel: The 0.1-200 keV count rate 
spectrum of XTE J1118+480, superposed with the 
best fit model consisting of an absorbed bb plus 
COMPTT in the spherical geometry (see text). Bot- 
tom panel: Ratio between the count spectrum and 
the model. 



perature of the seed photons, the best fit TVh is 
closer to -/Vjjg, even though the fit quality does 
not improve (x 2 /v = 279/227 or 280/227 ac- 
cording to the geometry of the electron cloud, a 
sphere or a slab, respectively; see Table 1 for the 
sphere). The cause of it is that the derived elec- 
tron temperature, kT = 239*gQ keV for a sphere 
and kT — 25ligg keV for a slab, is too much large 
to account for the observed cutoff at an energy of 
^100 keV (see also Figure 5). 

Then we tested the thermal- Compton model 
(compps v3.4 n ) of Poutancn & Svensson (1996). 
The advantage of this model is that, in addi- 
tion to allow for different geometries of the elec- 
tron cloud, it allows for the presence of Comp- 
ton reflection (Magdziarz & Zdziarski 1995) from 
a cold medium (like an accretion disk) with vari- 
able element abundances with respect to the so- 
lar ones. We considered spherical or slab ge- 
ometries with blackbody seed photons distributed 
homogeneously in the plasma and the Thomson 
optical depth, r, corresponding to the radius or 
the half-thickness, respectively. Initially we as- 
sumed reflection from a medium with solar abun- 
dances. For consistency, we also allowed for a nar- 
row Gaussian line (a = 0.01 keV) at 6.4 keV. In 
the spherical geometry and allowing the presence 



11 COMPPS is available on the internet at ftp://ftp.astro.su.se/ 
pub/juri/XSPEC/COMPPS. 
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Fig. 6. — Top panel: The 0.1-200 keV count rate 
spectrum of XTE J1118+480, superposed with the 
best fit model consisting of an absorbed bb plus 
COMPPS model in the spherical geometry and a 
solar metal abundance (see text). Bottom panel: 
Ratio between the count spectrum and the model. 

of additional blackbody photons (apart from those 
irradiating the plasma), we find x 2 l v — 255/225 
at kT = 58±g 2 keV, see Table 1. Thus the best 
fit provides a lower value of the electron tem- 
perature, that allows a better description of the 
observed high energy cutoff with respect to the 
COMPTT model 12 , see Figure 6. An F-test for 2 ad- 
ditional parameters provides a value of 4.6 x 10~ 4 
for the probability that the x 2 reduction from BB 
+ COMPTT to bb + COMPPS is due to chance. 

Similarly we cannot exclude the slab geome- 
try for the same model, that provides a statisti- 
cally similar fit, \ 2 l v = 263/225, even if the lat- 
ter model seems to show in correspondence of the 
150-200 keV data point a positive residual that is 
not observed in the case of the spherical geometry 
(see Figure 6). 

No evidence of relative reflection was found 
with the above model in both geometries, with a 
2<t upper limit of ~0.16 (see Table 1). This result, 
caused by the absence of an Fe K edge, contrasts 
other black black hole binaries in the hard state 
(e.g., Gierlinski et al. 1997). To check the robust- 
ness of the absence of reflection, we first consider 



12 Note that the values of kT obtained with COMPPS are much 
lower from that derived from COMPTT, which is caused by 
inaccuracies in the shape of the high-energy cutoff in the 
latter model, see, e.g., Zdziarski, Johnson, & Magdziarz 
(1996). 



the dependence on the inclination. We find no ev- 
idence of relative reflection even at i = 79° (about 
the maximum allowed by the data, §2.1) with a 
2cr upper limit of 0.26. Also no evidence of a nar- 
row Fe K line is found in the data, with a 2a up- 
per limit of the line intensity of 2.4 x 10~ 4 pho- 
tons cm~ 2 s _1 with energy from 6.4 to 6.9 keV, 
independently of the geometry adopted (see Ta- 
ble 1). 

A low metal abundance, Z, of the reflecting 
medium is expected for halo stars, like XTE 
J1118+480. Spectrally, it has an effect of reducing 
the K edge in the reflected spectral component 
without changing the shape of its high-energy 
cutoff. We find that allowing for Z < Z Q only 
slightly improves the fit bringing x 2 ' l v to 247/225 
(246/225) in the spherical (slab) geometry. Table 
1 shows the best fit parameters in both geometri- 
cal configurations. We are not capable to discrim- 
inate between these geometries. We now observe a 
statistically significant Compton reflection, simi- 
lar to the upper limits found when ZJZq was fixed 
to 1. We notice that the Comptonization parame- 
ters are slightly lower than those found with Cyg 
X-l in the hard state (Gierlinski et al. 1997), in 
the assumption of a single temperature isotropic 
model. The fit with variable Z is shown in Figure 
7, which also shows the remaining small residuals. 

We have also considered a multi-component 
model consisting of the COMPPS model in which 
the spectrum of seed photons peaks below the ob- 
served range (fcT see( j fixed at 10 eV), plus an ad- 
ditional bb component plus the narrow Gaussian 
profile (see above) at 6.4 keV. Such a low energy of 
seed photons is expected, e.g., when Comptoniza- 
tion of thermal synchrotron emission dominates 
(e.g. Wardzihski & Zdziarski 2000), see §6. This 
model, in both geometries (sphere or slab), yields 
a good description of the data as the last model, 
with a x 2 l v = 244/224. The best fit parameters 
are reported in Table 1. Notice the value of Nn 
in the slab geometry (= 1. 73^35 x 10 20 cnT 2 ) 
which is somewhat above -/Vh.g- 

Leaving the metallicity free to vary, we have 
also tested for the effect of replacing the black- 
body spectrum by a disk blackbody, in this case 
assuming the seed photons for Comptonization to 
be also those of the disk blackbody. This model 
yields plasma parameters very similar to the black- 
body case above and a similar fit quality (x 2 jv = 
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Fig. 7. — The BeppoSAX spectrum (crosses) fit- 
ted by (solid curve) a blackbody at 36 eV plus 
a thermal Comptonization in the slab geometry 
and a Compton reflection from a low-metallicity 
medium. The latter is shown separately by the 
dashed curve while the blackbody is shown by a 
dotted line. The spectrum is normalized to the 
level of the MECS (see text). The bottom panel 
shows the fit residuals. 
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Fig. 8. — The absorption-corrected model spec- 
trum (solid curve) corresponding to the fit shown 
in Fig. 7. The dashed curve shows the Comp- 
ton reflection component, while the dotted curve 
shows the blackbody spectrum observed. The ini- 
tial (i.e., before Compton scattering) spectrum of 
the seed photons incident on the plasma corre- 
sponds to the peak at ~0.15 keV. 



244/224 for a slab and = 245/224 for a sphere). 
The inner disk temperature is kT m = 41+4 e ^ 
at an inner radius greater than ^500 km (scaled 
at d = 1.5 kpc), that corresponds to ~ 34i? g 
(R g = GMx/c 2 ), assuming M x = 1OM . 

The absorbing columns fitted in models using 
COMPPS are close the Galactic TVh, with some ab- 
sorption intrinsic to the source possible as well. 
We note that the value of 2.8 x 10 20 cm~ 2 favored 
by Dubus et al. (2001) is inconsistent with our 
data, yielding, e.g., x 2 l v = 308/226 in the model 
with Comptonization in the spherical geometry 
This confirms our independent conclusions of §2.2. 

From the overall best-fit models, the mean flux 
level of the source during the observation is 5.7 x 
10~ 10 crg cm~ 2 s _1 , 6.9 x 10~ 10 erg cm~ 2 s _1 , 
1.9xl0 _9 erg cm~ 2 s _1 , and 3.4xl0~ 9 crg cm~ 2 s~ 
in the 0.1-2 keV, 2-10 keV, 15-200 keV, and 
0.1-200 keV bands, respectively (normalized to 
the MECS detector). The total (0.01-10 3 keV) 
flux corrected for absorption is somewhat model- 
dependent, and it equals to 4.2 xlO -9 erg cm~ 2 s _1 
for the BB + COMPPS model in the slab configu- 
ration (Z/Z sun free to vary in the fit) (see Table 
1 and fig. 7). The absorption-corrected flux of the 
observed blackbody is 2.2 x 10~ 10 erg cm' 2 s _1 to 
be compared with aflux of 3. 8xl0~ 10 erg cm~ 2 s _1 
due to blackbody photons incident on the hot 
plasma, i.e. 9% of the total flux (with an appar- 
ent amplification factor of 11). 

A temporal analysis of the data with Fourier 
techniques shows that the fractional variation of 
the source flux in the 0.01-50 Hz frequency range 
and the 2-10 keV energy band is 42% rms, which 
is consistent with the findings by Revnivtsev et 
al. (2000). In the 0.1-2 keV range, the fractional 
variation is even higher (62% rms). Also, a QPO 
with the centroid frequency at ~ 0.08 Hz is appar- 
ent in our power spectral density estimate. We do 
not find any statistically significant modulation at 
the optical period of 0.17 day. This implies that 
the X-ray source is axially symmetric and no ab- 
sorption in the line of sight is associated with the 
ouftlow of matter from the companion star. 

5. Spectral comparison with RXTE data 

RXTE performed many short pointings of XTE 
J1118+480. McClintock et al. (2001b) reported 
spectral results on this source based on the ob- 
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servation of April 18.8111-18.9583 UT and, for 
the study of the high energy part of the spectrum 
(>100 keV), results of the summed spectra ob- 
tained between April 13.39 UT and May 15.39 UT. 
No evidence of a high energy cut-off was found in 
the summed spectra, while the high energy data of 
the April 18 observation were affected by large un- 
certainties (see figure 4 of McClintock et al. 2001). 
One of the RXTE pointings of XTE J1118+480 
was performed simultaneously with a part of our 
observation. The RXTE observation (see, e.g., 
Revnivtsev et al. 2000) started 2000 April 15 
7:57:52 UT, and had an exposure time of 1152 s 
for the PCA, and 378 s and 361 s for the 2 High 
Energy X-ray Timing Experiment (HEXTE) clus- 
ters, respectively. A detailed study of this and 
other RXTE observations of XTE Jl 118+480 will 
be given elsewhere. We first fitted the RXTE 
spectrum by an e-foldcd power law, obtaining 
r ~ 1.8, i.e., a higher index than T ~ 1.7 for the 
BeppoSAX spectrum (Table 1). This difference 
is likely to be of instrumental origin, as a similar 
difference between the PCA and other X-ray in- 
struments is commonly found in fits to the Crab 
as well as other sources (e.g., Gicrlihski et al. 1999; 
Done, Madejski, & Zycki 2000). 

Then, we fit the 3-200 keV data with Comp- 
tonization in the spherical geometry, Compton re- 
flection at Zq, and a narrow Fe Ka line. We ob- 
tain J7/27T = 0.05 at the best fit, confirming the 
weakness of Compton reflection. However, the fit- 
ted electron temperature is somewhat larger than 
that of the BeppoSAX data, kT = 190±go k °V. 
We find, however, that if apply a correction of 
Ar = 0.1 to the model fitted to the PCA data, 
kT becomes equal to 140^ioo keV, consistent with 
the values of kT in Table 1 . The large confidence 
regions of kT are due to the limited statistics of the 
HEXTE data (the RXTE observation lasted only 
a small fraction of the BeppoSAX one). The fact 
that the summed RXTE spectra (covering more 
than a month elapsed time) are still described by 
a power law with no high energy cut-off below 
- 150 keV (McClintock et al. 2001) could be the 
result of the evolution of the source spectrum. In 
fact our preliminary results of the later 2 obser- 
vations of the source (see fig. 1) show a decreas- 
ing evidence of the high energy cut-off with time 
(Frontcraet al. 2001b). 

Summarizing, also the RXTE data appear to be 



compatible with the BeppoSAX ones, especially 
when the difference in the calibration between the 
PCA and BeppoSAX is taken into account. 

6. Theoretical Implications 

The Eddington luminosity at an H mass frac- 
tion of 0.7 can be written as 



L E ~ 1.5 x lO^Mx/lOM^ergs" 
and the gravitational radius is 



(8) 



R E = GM x /c 2 ~ 15(M X /1OM ) km. (9) 

The total X-ray luminosity corrected for absorp- 
tion calculated using the Comptonization model 
shown in Fig. 8 and assuming isotropy is 



1.1 x 10 36 ((i/1.5kpc) 2 ergs- 1 



8 x 10" 



1.5 kpc 



M x 
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(10) 



Here we used the estimates for the black hole 
mass and the source distance (§2.1). We note 
that the Lx above does not represent the bolo- 
metric luminosity of the source, whose EFe spec- 
trum has been observed by Hynes et al. (2000) to 
have a maximum in the far UV at a level similar 
to that of the peak at ~ 100 keV. Given their re- 
sults, the bolometric luminosity can be estimated 
as ~ 2Lx- Still, Lx represents the luminosity of 
the hot plasma flow. 

A fundamental issue is the geometry of the 
source. The weakness of Compton reflection 
and of blackbody emission poses some contraints. 
There are two possible solutions we can envis- 
age. The direct emission from the accretion disk 
is smeared by Compton scattering in the corona, 
or the hot gas subtends a small solid angle by a 
cold medium. 

A static disk/corona geometry (e.g., Haardt & 
Maraschi 1991), similary to the case of Cyg X-l 
in the hard state, can be reconciled with the data, 
as the large value of r/ cos i would scatter most of 
the reflection and blackbody emission in the power 
law continuum. However, because of the tight cou- 
pling of the hot and cold phases, this model cannot 
produce a spectrum as flat as that observed. In 
principle, a mildly relativistic coronal outflow can 
alleviate the problem (Beloborodov 1999; Malzac, 
Beloborodov, & Poutanen 2001). 
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Given that, the specific coronal model of Mer- 
loni, Di Matteo, & Fabian (2000) has some diffi- 
culties. It can produce a flat spectrum, but, as it 
postulates magnetic flares in a disk-corona geome- 
try, it predicts the presence of a strong blackbody 
component about an order of magnitude above 
that observed, even with relativistic bulk motion 
of the flares taken into account (compare their Fig. 
1 with our Fig. 8). 

On the other hand, the observed spectrum is 
fully compatible with Compton scattering in a hot 
central disk. Indeed, this is the geometry of our 
overall best-fit model. Then, the weakness of the 
blackbody and Compton-reflection components is 
naturally accounted for by a small solid angle sub- 
tended by any cold medium as seen from the hot 
disk. An important issue here is the origin of the 
seed photons. The weakness of the blackbody cor- 
responds to rather small size of the cold region. 
In either the blackbody or disk-blackbody model 
(see §4 and Table 1), the characteristic radius of 
the blackbody emitter is >^ 30i? g (at d = 1.5 
kpc, M = 10M Q ) for our best fit models. Thus, 
the seed photons can be supplied by an outer cold 
disk, probably overlapping with the hot one. 

Another compatible source of seed photons 
compatible is thermal synchrotron emission of the 
hot plasma itself (which naturally produces ho- 
mogenoeous and isotropic seed photons). Then 
the soft excess found in the data can represent 
either the peak of the self-absorbed thermal syn- 
chrotron emission (e.g. Wardzihski & Zdziarski 
2000) or blackbody emission not related to the 
dominant source of seed photons. The latter 
model was found in §4 to provide a good fit to the 
data. In fact, Wardzihski & Zdziarski (2000) have 
predicted that the thermal synchrotron process 
is able to dominate the supply of seed photons 
at an intermediate range of Eddington factors, 
Lx/Le < 10~ 3 , compatible with the case in XTE 
J1118+480. 

By considering the synchrotron process, we can 
constrain the magnetic field strength in the source. 
The turnover energy, E t , at which the thermal self- 
absorbed synchrotron emission peaks, is approxi- 
mately 



E t ~ 1.3 
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(Wardzihski & Zdziarski 2000). This yields a con- 



straint on the magnetic field strength of B < 10 8 
G, with the approximate equality corresponding to 
the observed soft excess at ~ 0.15 keV (Fig. 8) be- 
ing due to the turnover peak and no non-thermal 
electron tail. For lower values of B, the observed 
soft excess is due to blackbody emission. Then, 
the UV data of Hynes et al. (2000) constrain E t > 
10 eV. Wardzihski & Zdziarski (2001) have found 
that Et of equation (11) can increase by a factor 
of ~ 2-3 in the presence of a weak non-thermal 
tail beyond the Maxwellian electron distribution 
(with a corresponding high-energy photon tail ob- 
served in Cyg X-l, McConnell et al. 2000). Then, 
the UV constraint implies B > 10 7 G in the ab- 
sence of a non-thermal tail and B > 3 x 10 6 G 
with a tail. Only lower values of the above range 
of B correspond to equipartition in a hot accre- 
tion flow with moderately sub-virial ions [e.g., eq. 
(34) in Wardzihski & Zdziarski 2000], which indi- 
cates that the observed soft X-ray excess is due to 
a separate blackbody component rather than the 
turnover peak at E t . Summarizing, B < 10 8 G is 
implied directly by the data, and B > (3-10) x 10 6 
G for the sychrotron process to be the dominant 
source of seed photons. 

The high erratic time variability we find in the 
soft X-ray band (as well as in hard X-rays) points 
out to the origin of soft X-ray excess in the same 
spatial region as that occupied by the hot plasma. 
This is in agreement with both the seed pho- 
tons being provided by the blackbody emission of 
clouds within the hot plasma or their origin from 
the synchrotron process. The latter model also 
agrees with a similar time variability observed in 
the optical band, as noted by Merloni et al. (2000). 

Within the Compton scattering model in a hot 
central disk, our data appear compatible with 
the hot disk being advection-dominated, as re- 
cently proposed by Esin et al. (2001) on the basis 
of near-simultaneous HST, EUVE, Chandra and 
RXTE observations. However, if such interpreta- 
tion is correct, it is not clear how to compare this 
result with Cygnus X-l, which has a much larger 
luminosity and almost the same spectral parame- 
ters (Esin et al. 1998; Frontera et al. 2001a). 

McClintock et al. (2001b) reported a break 
at about 2 keV in the EF(E) spectrum of XTE 
J1118+480, which was interpreted by Esin et al. 
(2001) as due to the presence of a warm absorber. 
The data used by MCClintock et al. (2001b) were 
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obtained combining together different instruments 
aboard different missions. We do not see this 
break (see Figure 8), that could also be due to 
cross-calibration problems. 

7. Conclusions 

Our main findings can be summarized as fol- 
lows. 

We find the Galactic column density of Hi mea- 
sured in radio in the direction to the source as 
7V H ~ (1.281.44) x 10 20 cm" 2 (Dickey & Lock- 
man 1990) to be consistent with spectral fitting of 
our BeppoSAX data (see Table 1). On the other 
hand, we find the value of 2.8 x 10 20 cm -2 favored 
by Dubus et al. (2001) based on their measure- 
ment of Ca absorption to be intrinsically highly 
uncertain (§2.2) as well inconsistent with our data 
(§4). 

Our BeppoSAX spectrum of XTE J1118+480 
shows a power law with T ~ 1.7 modified by a soft 
X-ray excess and a high energy cutoff, with both 
features highly significant statistically. The soft 
excess contains a small fraction of the total flux, 
and is consistent with a blackbody component, ei- 
ther undergoing Comptonization or not. The high 
energy cutoff is well fitted by Comptonization in 
a thermal plasma with an electron temperature 
of ~ 10 2 kcV and a Thomson optical depth of 
t ~ 1. In agreement with the weakness of the 
blackbody, Compton reflection is weak, also con- 
sistent with the weakness of an Fe Ka line. The 
low amplitudes of all the features originating in 
the cold matter can have several explanations. If 
such components come from an underlying accre- 
tion disk, the relative large r and inclination angle 
easily dim the observable reflection and blackbody 
features by a factor exp(— r/cosi) ~ 0.1, consis- 
tent with our observations. However, in this case, 
the resulting Compton spectrum would be steeper 
than observed. A possibility is that any cold mat- 
ter subtends a small solid angle as seen by the hot 
gas. In this case the seed photons for Comptoniza- 
tion can be provided by an outer cold disk, or by 
synchrotron emission internal in the hot disk, or 
both. If the geometry is that of a central hot disk, 
it is not clear if our data are consistent with the 
ADAF solution, or instead if they rather indicate 
a radiative cooled hot accretion disk. 

The data are statistically consistent with the 



reflector to have a low metallicity, Z/Zq ~ 0.1, as 
expected due to the high latitute of the system. 
This could be an important finding, which should 
be tested by future observations of the secondary. 
The values of kT and r we find are very simi- 
lar to those typical of most luminous (e.g., Lx ~ 
0.02Le in Cyg X-l) black-hole binaries in the hard 
state, whereas the X-ray luminosity in our mea- 
surement is only ~ 10 _3 Le- 
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